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We use Monte Carlo simulations to investigate the ther- 
modynamical behaviour of aggregates consisting of few su- 
perparamagnetic particles in a colloidal suspension. The po- 
tential energy surface of this classical two-level system with 
a stable and a metastable "ring" and "chain" configuration 
is tunable by an external magnetic field and temperature. 
We determine the complex "phase diagram" of this intrigu- 
ing system and analyze thermodynamically the nature of the 
transition between the ring and the chain "phase" . 
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With progressing miniaturization of devices, there is 
growing interest in the thermodynamical behaviour of 
finite-size systems . A central question in this respect 
is, whether small systems can exhibit well-defined tran- 
sitions that could be interpreted as a signature of phase 
transitions which, strictly speaking, are well defined only 
in infinite systems. So far, reproducible features of the 
specific heat have been interpreted as indicators of "melt- 
ing" transitions in small rare gas clusters ||,|| . 

Here we investigate the intriguing thermodynamical 
behaviour of a structurally relaxed finite system which 
is controlled by two external variables, namely the tem- 
perature T and the magnetic field B cxt - The system 
of interest consists of few near-spherical, superparam- 
agnetic particles with a diameter of « 10 — 500 A in a 
colloidal suspension. Such systems, covered by a thin sur- 
factant layer, are readily available in macroscopic quan- 
tities, are called ferrofluids, and are known to form com- 
plex labyrinthine 0] or branched structures Q as many 
particle systems, whereas the only stable states for sys- 
tems with few particles (N < 14) are the "ring" and the 
"chain" §. 

The existence of two environmental variables, yet still 
only two isomer states, gives rise to a thermodynamical 
behaviour of unprecedented richness, as compared to that 
of other small clusters, such as the noble gas clusters 
[§,1). This is also an intriguing example of a classical, 
externally tunable finite two-level system. 

We will show that the system exhibits phase transitions 
between two ordered phases, one magnetic and the other 
one nonmagnetic, as well as phase transitions between 
these ordered phases and a disordered phase. Whereas 
the system is not susceptible to small magnetic fields, it 
shows a strong paramagnetic response when exposed to 
larger magnetic fields. 



Our model system consists of six spherical magnetite 
particles with a diameter of a = 200 A, mass m = 
1.31xl0 7 amu, and a large permanent magnetic moment 
fi = 1.68 xlO 5 he- The potential energy E p of this sys- 
tem in the external field B cxt consists of the interaction 
between each particle i and the applied field, given by 
Ui = —fti ■ -B cx t, and the pair- wise interaction between 
the particles i and j, given by || 

Uij = (MoAi ) [A* • Aj - 3(Ai ■ ryXAj ' %)] 



exp^-^— j-exp^-^— )j . (1) 

The first term in Eq. (Q) is the magnetic dipole-dipolc 
interaction energy. The second term describes a non- 
magnetic interaction between the surfactant covered tops 
in a ferrofluid that is repulsive at short range and at- 
tractive at long range We note that the most sig- 
nificant part of this interaction, which we describe by a 
Morse-type potential with parameters e = 0.121 eV and 
p = 2.5 A, is the short-range repulsion, since even at 
equilibrium distance the attractive part does not exceed 
10% of the dipole-dipole attraction. 

The equilibrium structures of small clusters are either 
rings or chains, both of which can be easily distinguished 
by their different mean magnetic moment (/z). 

In the following, we present the first complete "phase 
diagram" of a model aggregate of magnetic tops and de- 
scribe in detail how "phase transitions" occur in such a 
nanoscale system. Unlike in the bulk, where transitions 
between well-defined phases are sharp, small aggregates 
in a ferrofluid transform smoothly from one configuration 
to another due to changes in the two environmental vari- 
ables -Bext and T . Our results presented below are based 
on a careful analysis of set of 32 extensive Metropolis 
Monte Carlo simulations j7j , each of which consisting of 
6xl0 9 steps. 

The canonical partition function, from which all ther- 
modynamical quantities can be derived, is given by 



Z( J B ext ,T) = (27r/3)- 6A '/ 2 
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where (3 = (fceT) 1 and where the field B cxt is aligned 
with the z-axis. The pre-exponential factor addresses 
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FIG. 1. Monte Carlo results for the probability to find an aggregate in a state with its magnetic moment in the field 
direction [i z and potential energy E p . The individual contour plots show our results for the temperature T = 250 K at the 
field values (a) B cxt =0 G, (b) 73 cxt =40 G, (c) B cxt =60 G, and T =450 K at the field values (d) B cxt =0 G, (e) B cxt =40 G, (f) 
B oxt =60 G. 



the fact that each particle has three rotational and three 
center-of-mass degrees of freedom. The key quantities we 
monitor as a function of T and B cxt are the formation 
enthalpy E* of the isolated system that is decoupled from 
the external field, given by E* = J2%<j u ij = E p +(i z B cxt , 
and the z-component of the total magnetic moment of 
the aggregate, /j, z . By monitoring these two response 
quantities to the external temperature and field during all 
simulations, we determine the weighted density of states 

P{E* lf i z ;B cxtl T) = Z-^Sext.T) P (E* , p z ) 

x exp(-/?(£* - fi z B cxt )) , (3) 

which gives the probability to find the system in a state 
with the formation enthalpy E* and magnetic moment 
fi z . The partition function Z, which appears as the nor- 
malization constant, can be rewritten as 

Z(B ext ,T)= (2nf3)- eN / 2 (4) 
j dE*d(x z p{E*,fx z ) exp(-/3(£:* - ^ B cxt )) . 

We then combined the results of all simulations using the 
optimized Monte Carlo data analysis of Ferrenberg et al. 

in order to calculate the normalized density of states 
p{E*,fi z ) and p(E p ,fj, z ) with minimized statistical error 
|l0| ]. Using the above defined density of states, the field- 
and temperature dependence of the expectation value of 
any function F(E*,fj, z ) can be obtained as 

<F{E\p z -B cxt ,T) >= (5) 
J dfi z J dE* F(E* lf i z ) P{E*,fi z ;B ext ,T) . 

For the system described here we found it imperative to 
perform the simulations at sufficiently high temperatures 



in order to cover the whole configuration space properly. 
At low temperatures the thermal equilibrium might not 
be achieved even after extremely long iteration times, 
since transitions between rings and chains are very infre- 
quent and might never occur. In order to obtain a first 
idea about the stable and metastable states of the system, 
we plotted in Fig. [I] the probability P(E pi fx z \ B ext , T) of 
finding the aggregate in a state with potential energy E p 
and total magnetic moment in the field direction \i z . This 
is the projection of the probability to find the system in a 
specific state in the high-dimensional configuration space 
onto the (E p ,p z ) subspace. Regions in the (E p ,p z ) sub- 
space with a high probability indicate not only the ener- 
getic preference of corresponding states, but also a large 
associated phase space volume. 

Rings always have an absolute magnetic moment 
|^^max| tua ^. - g c i ose t zero. Consequently, also the z- 
component of the magnetic moment of rings is near zero, 
as seen in Fig. [j]. Even though the absolute magnetic mo- 
ment |/j,//j, max | of chains is close to one, these aggregates 
can not easily be distinguished from rings in the absence 
of a field. In zero field, chains have no orientational pref- 
erence and the z-component of their magnetic moment, 
/z z //x™ ax , also averages to zero. Of course, chains - un- 
like rings - do align with a nonzero magnetic field and, 
especially at low temperatures, show a magnetic moment 
/j, z //j,™ ax «l in the field direction. 

The relative stability of an aggregate is reflected in its 
potential energy E p . We find E p to increase (correspond- 
ing to decreasing stability) with increasing temperature. 
On the other hand, applying a magnetic field destabilizes 
rings in favour of field- aligned chains. With increasing 
field, chains are confined to a gradually decreasing frac- 
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tion of the configurational space which sharpens their 
distribution in the (£p, subspace, as seen when com- 
paring Figs. |](a)-(c) and Figs. |](d)-(f). 

Under all conditions, we find two more or less pro- 
nounced local maxima in the probability distribution P, 
corresponding to a ring with 0</i z //z™ ax <i;l, and a chain 
with 0^/i z //i™ ax <I. At zero field we observe a predom- 
inant occupation of the more stable ring state. Due to 
the relatively small energy difference with respect to the 
less favourable chain AE" /N = (E^ hain - E^)/N = 
0.06 eV, both states get more evenly occupied at higher 
temperatures. At fields as low as B cxt = 40 G, the energy 
difference between chains and rings drops significantly to 
AE^/N = 0.02 eV. As seen in Fig. |(b), this results in 
an equal occupation of both states even at low tempera- 
tures. At the much higher field value i? ox t = 60 G, chains 
are favoured with respect to the rings by a considerable 
amount of energy AE" /N = —0.2 eV. This strongly sup- 
presses the occurrence of rings, as seen in Figs. |](c) and 

(f)- _ 

Fig. [j] shows not only the stable and metastable states 
under given conditions, but also the states found along 
the preferential transition pathway between a ring and 
a chain in the projected (E p ,fx z ) subspace. During the 
transition between a chain and a ring, each aggregate 
must undergo a continuous change of E p and [i z . The 
favoured transition pathways are then associated with 
high-probability trajectories in the (E p ,/j, z ) subspace. 
The value of the activation barrier AE act is then given by 
the largest increase of E p along the optimum transition 
path which connects the stable and metastable ring and 
chain islands. In our simulations we found that the ac- 
tivation barrier occurred always at /U z //Lt™ aa: «0.22. Con- 
sequently, we concluded that the field dependence of the 
activation energy follows the expression A_E act (i? cx t) = 
AE; ct (B cxt = 0) - 0.22 /^ ax S ext . 

In order to quantitatively describe the "phase transi- 
tions" occurring in this system, we focused our attention 
on the specific heat and the magnetic susceptibility. The 
specific heat per particle in a canonical ensemble is given 
by cb — d(E/N)/dT, where the total energy is given by 
E = ^Nk B T + E p . Correspondingly, we define the mag- 
netic susceptibility per particle as \ = d(/x z /iV)/d_B cxt . 
These response functions are related to the fluctuations 
of -Ep and \i z by 



cb 



6N 



k B + k B /3 2 ((E 2 ) - (E) 2 ) 



/N, 



X =[(3((£)-(» z ) 2 )]/N 



(6) 
(7) 



Phase transitions are only well-defined in infinite systems 
and are associated with a discontinuous change in the 
total energy and specific heat when crossing the phase 
boundary. The corresponding changes in finite systems 
are more gradual. This is seen in the well defined, yet 
not sharp "crest line" separating the ring and the chain 
"phase" in the T — B cxt "phase diagram" in Fig. ||(a). 
These results illustrate how the critical magnetic field for 
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FIG. 2. Specific heat per particle cb of the system as 
a function of temperature T and the external magnetic field 
Bcxt- Results for the entire temperature and field range in- 
vestigated here are presented as a contour plot in (a). The 
temperature dependence of cb for selected values of B e xt is 
presented in (b). 



the ring-chain transition decreases with increasing tem- 
perature. At high temperatures, the "line" separating 
the "phases" broadens significantly into a region where 
rings and chains coexist. 

The line plot in Fig. ||(b) is the respective constant- 
field cut through the contour plot in Fig. ||(a). As can 
be seen in Fig. ||(b), there is no transition from chains 
to rings, indicated by a peak in cb, at fields exceeding 
50 G, which is close to the critical field value at which 
chains become favoured over rings at zero temperature. 
At fields B cxt << 40 G, on the other hand, there is no re- 
gion where chains would be thermodynamically preferred 
over the rings, and we only observe a gradual transi- 
tion from the ring phase into the coexistence region with 
increasing temperature. The specific heat behaviour at 
zero field resembles that of a small system with a gradual 
melting transition close to 150 K and an onset of disor- 
der at about 350 K p]| . As seen in Fig. ||(b), the critical 
temperature and the width of the transition region can 
be externally tuned by the second thermodynamical vari- 
able, the external magnetic field B cxt - 

Fig. H displays the magnetic susceptibility x> another 
prominent indicator of phase transitions in magnetic sys- 
tems, as a function of T and i? ox t. Like the specific 
heat in Fig. ||(a), the crest line in \ separates the chain 
"phase" from the ring "phase" in this T — B ext "phase di- 
agram" . Moreover, Fig. |3| reveals the fundamentally dif- 
ferent magnetic character of these "phases" . Whereas the 
system is nonmagnetic in the ring "phase" found below 
40 G, it behaves like a ferromagnet consisting of Langevin 
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FIG. 3. Magnetic susceptibility per particle \ °f the sys- 
tem as a function of temperature T and the external magnetic 
field Bcxt- Results for the entire temperature and field range 
investigated here are presented as a contour plot in (a). The 
temperature dependence of x f° r selected values of B eyL t is 
presented in (b). 



paramagnets in the chain "phase" at higher fields. The 
transition between these states is again gradual, as ex- 
pected for finite systems. The line plot in Fig. ||(b) is the 
respective constant-field cut through the contour plot in 
Fig. ||(a). When the system is in the chain "phase" it 
behaves like a paramagnet obeying the Curie- Weiss law, 
as can be seen in Fig. ||(b) fl^| . 

At relatively low temperatures, where the aggregates 
are intact, the expectation value of the magnetic mo- 
ment first increases due to the gradual conversion from 
nonmagnetic rings to paramagnetic chains. According 
to Fig. |(b), this uncommon behaviour persists up to 
T = 200 K at B cxt = 40 G. This trend is reversed 
at higher temperatures, where all aggregates eventually 
fragment into single paramagnetic tops. In this temper- 
ature range, the magnetic moment as well as the suscep- 
tibility decreases with increasing temperature. 

Since the transition probability between both states is 
extremely low at low temperatures and fields, magneti- 
cally distinguishable metastable states can be frozen in. 
A chain configuration, which is metastable in zero field, 
can be prepared by first annealing the system to T>350 K 
and subsequent quenching in a strong field. Similarly, a 
frozen-in ring configuration is unlikely to transform to a 
chain at low temperatures, unless exposed to very large 
fields. Thus the above described phase diagrams can be 
used to externally manipulate the self-assembly of mag- 
netic nanostructures. 

In conclusion we have studied the thermodynamical 
behaviour of a finite two-level system, which is externally 



tunable by two independent variables, namely the tem- 
perature and the magnetic field. Much of the behaviour 
encountered in this system, such as transitions between 
different states, has a well-defined counterpart in infinite 
systems. The reason for the encountered richness of the 
thermodynamic and magnetic properties is the relative 
ease of structural transformations, which is typical for fi- 
nite systems. Consequently, we expect other finite mag- 
netic systems, e.g. small transition metal clusters, where 
a small number of structural isomers with substantially 
different magnetic moments could coexist |l3[ ], to follow 
this behaviour. Moreover, we expect that our results 
can also be transferred to nanocrystalline material, such 
as magnetic clusters encapsulated in the supercages of 
zeolites, which will likely retain some of the intriguing 
properties of the isolated finite systems. 
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